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I ^1 Abstract. The lattice dynamics of coesite has been studied by a combination of 

^ diffuse x-ray scattering, inelastic x-ray scattering and an ab initio lattice dynamics 

^ calculation. The combined technique gives access to the full lattice dynamics in 

harmonic description and thus eventually provides detailed information on the elastic 

^^ properties, the stability and metastability of crystalline systems. The experimentally 
CO 

' validated calculation was used for the investigation of eigenvectors, mode character 
■^ 

(^ and their influence on the density of vibrational states. High symmetry sections of 

. the reciprocal space distribution of diffuse scattering and inelastic x-ray scattering 

^^ spectra as well as the density of vibrational states and the dispersion relation are 

reported and compared to the calculation. A critical point at the zone boundary is 



^ 
^ 



X 



Vh found to contribute strongly to the main peak of the low energy part in the density of 

vibrational states. Comparison with the most abundant Si02 polymorph - a-quartz - 
reveals similarities and distinct differences in the low-energy vibrational properties. 



PACS numbers: 63.20.D-, 63.20.dd, 63.20.dk, 72.10.Di 
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1. Introduction 

Coesite, Si02, is the highest density tetrahedrally coordinated crystalhne polymorph of 
sihca with space group C 1 2/c 1 [Ij. It was first synthesized in 1953 at 3.5 GPa and 
750 °C [2] and later found in sandstone of the Arizona Barringer crater |3j and led to 
the general acceptance of the impact cratering theory and to important implications 
for the recognition of meteorite impact craters in quartz-bearing geologic formations 
[1] . The discovery of coesite revolutionized the whole meteorite study and the analysis 
of impact products. It was also shown that being encapsulated in diamond formed 
in deep Earth's interior, coesite allows the unambiguous identification of 'fossilized' 
high pressure states in individual inclusions of mantle samples delivered to the Earth's 
surface in the course of various geological processes [5j| . The fossilized pressure phases of 
the inclusion and the thermoelasticity of the host-inclusion (diamond-coesite) ensemble 
may provide a highly accurate geobarometer [5]. Elastic and dynamical properties of 
coesite have been intensively studied in the past. Elastic constants for instance were 
examined using Brillouin light scattering [6] and ab initio calculations [7]; vibrational 
properties at zero momentum transfer were measured by Raman [HI E] and infrared 
spectroscopy fiUl \TT\ . The compression mechanism was studied by single crystal x- 
ray diffraction [121 [13] and ab initio calculations [7] . The investigation of silica under 
pressures lead to the discovery of unexpected phenomena, such as, for example, the 
formation of mesoporous coesite at a pressure of 12 GPa and a temperature of 300 °C 
[m [15]. Currently, the study of the phonon dispersion relations is limited to numerical 
calculations using potential based methods [16j. Experimental studies of dispersion 
relations have not been presented yet. The knowledge of the full lattice dynamics is, 
however, fundamental for the understanding of the compression mechanism and phase 
transitions. In particular the origin of the first peak of the density of vibrational states 
(VDOS) is of interest regarding the origin of the Boson peak in silica glass. As suggested 
in Ref. [TTj the Boson peak in glasses originates from the acoustic phonon branches 
near the boundary of the pseudo-Brillouin-zone and has its counterpart in the VDOS 
of the corresponding crystal. In the following we report the results of the powerful 
combination of diffuse scattering and inelastic x-ray scattering (IXS) for the study of 
distinct lattice dynamical features in a twinned crystal and the validation of the DFPT 
calculation. The combined approach allows understanding the lattice dynamics at 
arbitrary momentum transfers over the entire energy range and thus delivers a complete 
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picture. The validated calculation is used for the study of particularities in the dispersion 
relations and contributions of different modes to the VDOS. The paper is structured as 
follows. The synthesis of the high quality crystal and the details of the measurement 
are given in Experimental methods. The lattice dynamics calculation is described in 
section Calculation. Experimental and theoretical results are presented and discussed 
in Results and Discussion. The Conclusion summarizes our findings. 

2. Experimental methods 

A polycrystalline sample of coesite was synthesized at the university of Bayreuth using 
the high-pressure high-temperature technique at P = 5.5 GPa and T = 1000 °C. A 1000- 
ton hydraulic press (Voggenreiter GmbH) with a toroidal-type high-pressure cell |18] was 
employed. The size of the synthesized sample was of about 50-70 mm^ of pure coesite. 
As a starting material compacted amorphous silica (Sigma Aldrich, 99.9 % purity) 
was used. The cross-section of the sample container just after the synthesis is shown 
in Figure [l] The single crystals were grown by the hydrothermal method described 
elsewhere |19]. The x-ray diffuse scattering experiment was conducted on beamline 
ID29 [20] at the European Synchrotron Radiation Facility (ESRF). Monochromatic X- 
rays with wavelength 0.7 A were scattered from an elongated (^ 0.1 x 0.3 mm^) coesite 
crystal at room temperature in transmission geometry. The sample was rotated with 
an increment of 0.1° orthogonal to the beam direction over an angular range of 360° 
while diffuse scattering patterns were recorded in shutterless mode with a PILATUS 
6M detector [21j . The orientation matrix and geometry of the experiment were refined 
using the CrysAlis [22] software package. The sample was found to be twinned in the 
ac plane. The single crystal IXS study was carried out on beamline ID28 at the ESRF. 
The spectrometer was operated at 17.794 keV incident energy, providing an energy 
resolution of 3.0 meV full-width-half-maximum with a beam focus of 30 x 60 /im. IXS 
scans were performed in transmission geometry along selected directions in reciprocal 
space. Further details of the experimental set-up and the data treatment can be found 
elsewhere [23] . 

The generalized x-ray weighted VDOS (X-VDOS) was obtained from IXS spectra 
of a polycrystalline sample measured at ID28 [24j. The scattered radiation was collected 
by nine crystal analysers. The momentum transfer resolution of each analyser was ~ 
0.3 nm~^. The values of the momentum transfers for each analyser were chosen away 
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from the Debye-Scherrer rings and covered the [10 : 70] nm~^ range. The data combine 
the results of measurements with 1.4 meV resolution at 23.725 keV incident energy 
within [-25 : +25] meV energy range and 0.2 meV energy steps and results from the 
measurement with 3.0 meV resolution within [-25 : +180] meV energy range and 0.7 
meV steps. The elastic peak in the IXS spectra was subtracted using the instrumental 
function of each analyser determined by x-ray scattering from a polymethylmethacrylate 
(PMMA) sample close to the maximum of its structure factor. The X-VDOS was 
obtained from the summed IXS spectra within the incoherent approximation following 
the data treatment procedure established in [21]. 

3. Calculation 

The lattice dynamics calculation was performed using the DFPT approach |25] 
as implemented in the CASTEP code [26l |27] . The local density approximation 
within the plane-wave pseudopotential formalism was employed using norm conserving 
pseudopotentials. The plane wave cut-off was set to 800 eV and the electronic 
structure was computed using a 3 x 3 x 2 Monkhorst-Pack grid. The structure 
was optimized with the Broyden-Fletcher-Goldfarb-Shannon method, which is a full 
geometry optimisation of lattice and internal parameters. The cell parameters of the 
optimized cell are compared to the experimental values [2B] from single-crystal neutron 
and x-ray diffraction in Table [ij They agree within 1.6 %. 

Phonon frequencies and eigenvectors were computed on a 4 x 4 x 3 Monkhorst- 
Pack grid of the irreducible part of the Brillouin zone by a perturbation calculation 
and further Fourier interpolated for the VDOS and dispersion relations. The calculated 
phonon energies were tested to be converged to < 0.05 meV. Thermal diffuse scattering 
(TDS) and IXS intensities were calculated from the phonon eigenvectors and frequencies 
following the previously established formalism p3| [221 ED] , assuming the validity of both 
the harmonic and adiabatic approximation. TDS intensities were calculated in first order 
approximation. 

4. Results and discussion 

High symmetry reciprocal space sections of diffuse scattering as obtained from 
experiment and corresponding calculated TDS intensity distributions are shown in 
Figure [2J Corrections for polarisation and projection [25 and the Laue symmetry 
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of the system were applied. A complex intensity distribution is noticeable. Intense 
features indicated in Figure |2] were selected for a detailed IXS measurement in order 
to distinguish between possible elastic and inelastic contributions and to resolve the 
energies of the phonons contributing most to the TDS. A remarkable consistence between 
the experimental and calculated patterns can be seen for instance from the shape of 
indicated features. The twinning is mostly visible in the third column in Figure |2} The 
number of Bragg reflections can only be described by a combination of two domains. 
The HKO plane is common. The structure could be solved for a twinning of 180° around 
the reciprocal lattice axis a* with an intensity contribution of 25% of the smaller crystal. 
The corresponding calculated TDS maps in the second row was created by the weighted 
superposition of the indicated intensity distributions. Experimental artefacts due to 
non-uniform absorption of the sample are visible, arising from the anisotropic shape of 
the sample. 

Figure [3] shows two representative IXS spectra. One can clearly note the influence 
of the two crystalline domains on the inelastic spectra in Fig. [sla). The spectra of 
the two domains belong to different Q-values with different phonon eigenvectors and 
energies, and are thus different. The calculated spectra of the two domains in the 
common plane Is] b) are, except for the overall intensity, identical. The experimental 
spectra also show a small elastic line. After scaling the calculated energies by 1.045 the 
theoretical spectra reproduce quite well both position and intensity of the phonons. The 
scaling factor was determined from the VDOS, and its value is justified further below. 
Figure |4] shows the IXS intensity maps along the indicated directions from calculation 
and measurement. Each experimental map consists of 4 measured spectra with linear Q- 
spacing and energy steps of 0.7 meV. The momentum- and energy- transfers are linearly 
interpolated to 20 q-points and 72 energy steps. The theoretical maps were obtained 
by the weighted superposition of inelastic spectra for the two domains. The inelastic 
intensity is calculated from the eigenvectors and eigenfrequencies for 120 points along 
the given direction in reciprocal space and convoluted with the experimental resolution 
function of the spectrometer. The experimental spectra show that the diffuse scattering 
is of almost exclusively inelastic nature. Taking into account the intensity contribution of 
the two crystalline domains, both, energy transfer and inelastic intensities of the different 
phonon branches are well reproduced by the calculation for arbitrary directions. This 
implies that the theory correctly predicts both eigenvalues and eigenvectors at arbitrary 
momentum transfers. 
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The generalized x-ray weighted VDOS (X-VDOS), 

G{E) = Ud'Qj2\Yl ^e^^^-'^-lQ • ec,,r.mE - Eq 



with momentum transfer Q, phonon energy E and eigenvector component e of 
branch j, atomic form factor /, mass m, position vector r, Debye Waller factor W of 
atom n and normalization factor A^ = Xln fn{Q) / y/frtn for coesite is shown in Figure 
5] a). The X-VDOS was used for the determination of an overall energy scaling factor. 
In fact the VDOS probes the energy of the ensemble of states in three dimensional 
reciprocal space and is therefore most significant for the determination of the scaling 
factor. It turns out that an overall stretching of 1.045 leads to a good agreement of 
experimental and theoretical X-VDOS. Discrepancies between experiment and theory 
are mostly due to the limited accuracy of sampling the reciprocal space with powder IXS 
spectra. The underestimation of the calculated energies can be attributed to the limited 
accuracy of the exchange correlation function within the local density approximation; 
see Ref. P?] for a detailed discussion. The partial density of states (Figure [SJb and c)) 
separate the contribution of silicon and oxygen atoms. Looking at the low energy part 
of the partial VDOS (Figure Isle)) we find that the first peak located at 10.1 meV has 
equal contributions from silicon and oxygen atoms. The main peak of the scattering 
function at 14.4 meV is, however, dominated by the vibration of the oxygen atoms. The 
low energy VDOS of a-quartz is plotted in Figure [s] f ) for comparison. The a-quartz 
calculation was performed using the same pseudopotentials and similar parameters as 
for coesite. For a detailed discussion on the calculation see Ref. [27]. Here the first 
peak dominates the low energy part of the VDOS. It arises as in coesite from equal 
contributions from silicon and oxygen atoms. 

Peaks in the VDOS require critical points in the phonon dispersion relations for 
which the gradient in all crystallographic directions becomes zero [32j. In order to 
localize the critical points contributing most to the first and the main peak in the low 
energy part of the VDOS, two filters were applied simultaneously. An energy filter of 
A(jj= 0.3 meV was applied to the ah initio calculated phonon eigenfrequencies of the 
first Brillouin zone and l/|VqiX'(q)| was computed within this energy window. A saddle 
point close to the zone boundary at (0.87 0.69 0.42) was found to contribute strongly 
to the first peak in the VDOS and a local minimum at Y (1 0) is responsible for 
the main peak of the VDOS. The local contribution of the critical points to the VDOS 
within Qc ± O.lr./.M. is shown in Fig. [sld). The dispersion relations along the main 
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crystallographic direction through the critical points and the displacement patterns are 
shown in Figure |6] and compared to the ones of a-quartz. 

Summarizing the results shown in Figure p\ and p\ we note: (i) Critical points close 
to or at the zone boundary contribute most to the first and main peak of the low energy 
VDOS in both systems. This observation might be explained in a simplified picture 
with the piling up of vibrational states due to a flattening of the dispersion relations at 
the zone boundary p2]. (ii) The first peaks in both coesite and in a-quartz are due to 
an almost equal contribution from silicon and oxygen atoms. The main peak in coesite 
is, however, dominated by the vibration of oxygen atoms. The character of vibration 
is more libration like, (iii) In coesite both the first and the main peak are located at 
higher energies than in a-quartz. This observation could be expected from the higher 
density structure, but we observe that (iv) the atomic displacements of both critical 
points in coesite are different from the atomic displacement pattern of a-quartz. This 
shows that the peaks cannot be compared directly. 

Calculated dispersion relations along high symmetry directions and the result of an 
interatomic potential calculation [16j together with experimental values from the IXS 
measurements and Raman scattering p] are shown in Figure^ The two calculations are 
in reasonable qualitative agreement for most of the branches. For the Y - F - direction 
our calculation agrees much better with the experimental results, particularly for the 
fast acoustic and the measured optic mode. The calculation of the low energy optical 
branches is very delicate and sensitive to small geometrical differences and require a 
fine electronic grid sampling. We note that both acoustic and optical phonon branches 
contribute to the main peak of the scattering function at 14.4 meV. The M - Y dispersion 
relation of the branch containing the critical point is very fiat. Experimental phonon 
energies as determined by the IXS measurement and Raman scattering [8] are in good 
agreement with our calculation. 

The calculation is further compared to experimental results from Raman and 
infrared [T^ studies in Table [2J The modes were attributed by a careful comparison 
of experimental and calculated intensities. A good overall accordance is obtained. The 
calculation predicts some additional modes to be Raman or infrared active with low 
intensity contribution. 
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5. Conclusions 

With tlie validated lattice dynamics calculation it is possible to access the dynamical 
properties at arbitrary momentum transfer. The investigation of the VDOS, probing the 
ensemble of vibrational states, allows the identification of distinct dynamical features. 
In the case of coesite a small linear energy scaling of the ab initio calculated phonon 
frequencies leads to a good agreement of experiment and theory. The investigation of 
the nature of the dominating features in the low energy part of the VDOS shows that 
critical points located close to or at the zone boundary provide the largest contribution. 
The contribution of oxygen atoms is found to dominate the main peak. The extension 
of the model for the lattice dynamics at ambient conditions to high pressure potentially 
builds the basis for understanding the compression mechanism and phase stability. This 
extension may provide an accurate model for the elastic properties at the formation of 
coesite inclusions in diamond implying a precise calibration of the coesite-in-diamond 
barometer proposed in Ref. ^. Comparison with the most abundant silica polymorph 
a-quartz reveals valuable new insights into the low-energy vibrational properties of this 
prototypical oxide. Within the framework of silica polymorphs we are now able to 
provide accurate models of the lattice djTiamics of a-quartz [30j, coesite and stishovite 
|33j . The calculation for these polymorphs can be extended to high pressures and allow 
the derivation of thermodynamical properties. The prediction of the lattice dynamics 
for new high pressure phases like seifertite and others |34j should be possible within the 
employed calculation scheme. 
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Table 1. Cell parameters of Si02 coesite. 
Calculation Experiment [28] 



a= 7.137 A 
b = 12.295 A 
c = 7.072 A 
a = 7 = 90° 
/3 = 120.374° 



a ■■ 
b 

C : 

a 



7.136 A 
12.384 A 

7.186 A 
7 = 90° 
120.375° 




Figure 1. (colour online) Cross-section of the pressure pellet containing a 

polycrystalline sample of coesite synthesised at P= 5.5 GPa and T=1000 C. 
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Figure 2. Experimental (first row) and calculated (second row) diffuse scattering 
intensity distributions of coesite in high-symmetry reciprocal space sections. The 
absolute intensity is scaled to the best visualisation of the diffuse features. The 
calculated intensity distributions in the second row are created by the weighted 
superposition of the two crystalline orientations shown in the third and fourth row. 
The diffuse features marked by circle and polygon were selected for an IXS study, the 
intense features highlighted by rectangle and polygon are guides for the comparison of 
experimental and calculated diffuse scattering. See text for further details. 
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Figure 3. (colour online) Experimental (black points with errorbars) and theoretical 
(red line) IXS spectra from coesite crystal. The theoretical spectra are the weighted 
sum of the two crystalline orientations: domain 1 - 75% (blue dashed line), domain 
2 - 25% (green dotted dashed line). Theoretical intensities were convoluted with the 
experimental resolution function and the energy transfer was scaled by 1.045. 
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Figure 4. Experimental IXS intensity maps (first row) from coesite crystal together 
with theoretical intensity maps (second row) along the indicated directions. See text 
for further details. 
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Figure 5. (colour online) a) Experimental (blue points) and theoretical (black line) 
X-VDOS of coesite. The calculated X-VDOS is the x-ray weighted sum of the partial 
VDOS of oxygen (b) and silicon (c). The energies of the theoretical data were scaled 
by 1.045 and the resulting spectrum was broadened by the experimental resolution, 
d) Low-energy part of the VDOS (black line) and the local contribution of the critical 
points to the VDOS within qc ± O.lr.l.u. (purple line). The partial VDOS of oxygen 
(red dots) and silicon (green line) are compared in (e), where the silica contribution is 
multiplied by a factor of two. f) The partial VDOS of oxygen (red dots) and silicon 
(green line) in the low-energy range of a-quartz. The silica contribution is multiplied 
by a factor of two. 
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Figure 6. (colour online) a) Dispersion relations along the main crystallographic 
directions through the critical point (0.87 0.69 0.42) and the displacement pattern 
(primitive cell) of the atoms at this saddle point contributing strongly to the first peak 
in the VDOS of coesite. b) Dispersion relations and displacement pattern at the Y 
point (1 0) in coesite. c) Orthogonal dispersion relations and displacement pattern 
at the L point (1/2 1/2) in a-quartz. 
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Figure 7. (colour online) Calculated dispersion relations along high symmetry 
directions of coesite crystal (black lines) and the results of an interatomic potential 
calculation [T5] (dashed green lines) together with experimental values from the IXS 
measurements (blue squares) and Raman measurements [H] (red points). The critical 
point at Y is marked with a magenta star. The energies of our calculation are scaled 
by 1.045 whereas the energies of interatomic potential calculation are scaled arbitrary 
for best visual fit. 
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